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Abstract 

The ground state Lamb shift of a semiconductor spherical Quantum Dot is computed in the effective mass 
approximation. It appears to be significant enough to be detectable for a wide range of small Quantum Dots 
synthesized in semiconductors. A possible way to observe it, via the Casimir effect, is suggested. 
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1. Introduction 

The investigation of atomic physics properties of semiconducting Quantum Dots (QDs) is a popular 
topic for fundamental as well as applied physics. Such nanostructures restrict the motion of charge carriers 
to a confined region of space. Two QDs are never identical, because of the crucial role of phonons, surface 
effects or bulk disorder on their electronic structure, but they may be advantageously considered as artificial 
giant atoms. Thanks to their adjustable quantized energy spectrum, controlled only by their size, they 
may be used, for example, as semiconductor lasers [l|, single-photon sources 0, qubits 0, single-electron 
transistors artificial fluorophores In the early 1980s, the so-called Quantum Size Effects (QSE), 
characterized by a blue-shift of their optical spectra, has been observed in a large range of strongly confined 
systems It comes from a widening of the semiconductor optical band gap, due to the increase 



of the charge carriers confinement energy [111 ]. Modern methods as well as empirical pseudo-potential 
methods or ab initio approaches, such as the Density Functional Theory (DFT), which are appropriate for 
numerical determination of semiconductor bulk band structures or confinement energies, arc discussed in 
Refs. [Sill. 

However, to the best of our knowledge, no comprehensive approach, which offers a significant 
analytic treatment, seems to exist at present. 

In this work, to explore their analogy with real atoms, we look at the interaction of a spherical semi- 
conductor QD with a quantized electromagnetic field and investigate the Lamb shift of its energy levels. 



Discovered and actively studied in the late 1940s [14| - |17|. th e Lamb shift still remains nowadays an area of 



intense research, for example within dielectric materials [18j , on the excitation modes of an electromagnetic 
cavity [lH, or on the coupling of electromagnetic modes to the semiconductor QD surface j2£j. Recent 
works also present experimental protocols, in which Lamb shifts enhanced by an electromagnetic cavity are 
measured in nanostructures like transmon [2l[ or studied in semiconducting QDs coupled to a planar slab 



of negative-index material |22j. However, the Lamb shift, well understood for real atoms, seems not to be 
approached neither theoretically nor experimentally in semiconductor QDs. 

The object of this paper is to fill this gap by offering a consistent theoretical treatment leading to the 
evaluation of the shifts for a wide range of realistic spherical semiconductor QDs. To describe a confined 
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electron-hole pair, it seems appropriate to use the popular effective mass approximation (EMA) model [231 — 
[2IH . We first review the tools needed to work out the Lamb shift in section [2j and the basic properties 
of a spherical semiconductor QD in section [3] From there on, we deduce, in sections [4] and [5j the Lamb 
shift for a massive charged particle confined by an infinite potential well. As this Lamb shift is negative, we 
comment, in section [6j our result to emphasize the difference with the positive Lamb shift observed in real 
atoms. In section [7J for spherical semiconductor QDs, we show that there exists a possibility to observe 
experimentally the Lamb shift, at least in a so-called strong confinement regime and for a judiciously chosen 
semiconductor. We also suggest a protocol, built on the Casimir effect, which might allow to observe this 
predicted Lamb shift, despite the fundamental non-degeneracy of semiconducting QD levels. A concluding 
section summarizes our main results and indicates possible future research perspectives. 



2. Considerations on Lamb shift 

In 1947, investigations on Zeeman effect in hydrogen atom showed that its fine structure docs not agree 



with the predictions of Dirac theory [14J. The 2s-level, which should be degenerated with 2p- levels, is 



actually shifted by an energy ps 1057MHz, the so-called Lamb shift. Theoretical work demonstrated that it 



arises from the coupling of the electron motion with the surrounding quantized electromagnetic field 15|, [16[ . 
Consider a non-relativistic spinless massive particle of mass m* and of charge qe interacting with a quantized 
electromagnetic fielcQ- It is represented by the Pauli-Fierz Hamiltonian in the Coulomb gauge [2^|, here 
written in units where h = c = 1, 

Hpp = i/ + ei/i n t + e 2 H' int + H em . 

2 

The particle Hamiltonian in absence of electromagnetic field Hq = -^r+V(r) is supposed to be diagonal in an 
orthonormal basis of its eigenvectors {|n)}„ with energy eigenvalues {-E^ln- Hint = — Q^?- an d H[ nt = q 2 ^-r 
are interaction Hamiltonians, and H cm = J d 3 r E + B is the free electromagnetic field Hamiltonian, A being 
the vector potential, E and B the electric and magnetic fields. The electromagnetic field is second quantized 
in the Coulomb gauge [3(|] . The interaction Hamiltonian H[ nt will not be taken into account on the basis 
of weak intensities of light sources involved^. We now review the results of the two main approaches to the 
Lamb shift. 

The Bethe approach to the Lamb effect is purely perturbative. The quantum second order time inde- 
pendent degenerate perturbation theory is applied to the Pauli-Fierz Hamiltonian Hpp , where Hamiltonians 
H mt and H[ nt are treated as perturbations in the weak field limit flBT ]. Using renormalization arguments, 
the Lamb shift for an energy eigenstate |n) is found to be given by 

^ = t&^ii^m) ? i<-w>i"«- - m 



-3,,„. i 1 H(E^i)J <n|vv(r)l " > ' (1) 
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where a = |- is the fine structure constant written in the chosen units system. (\E m — E n \) is the mean 
value of energy level difference absolute values. Historically, this predicts a Lamb shift for the hydrogen 



atom 2s-lcvel in excellent agreement with experimental values 15l. |17|. 

In the Wclton approach, the Lamb shift is interpreted as a fluctuation effect on the particle position due 
to its interaction with the surrounding electromagnetic field. These fluctuations Ar can be described as a 
continuous random variable, whose probability density is a three-dimensional centered isotropic Gaussian 
distribution of variance 

<(Arf> = -^log- V (2) 



1 Rclativistic corrections can generally be calculated [27I. l28l . but they are not relevant in the comprehension of Lamb effect, 
particularly in spherical semiconductor QDs, where confined charged carriers are considered as non-rclativistic particles. 
2 For more details, one can refer to complement Axm of chapter XIII in Ref. l3ll . 
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where k is a IR cut-off, and m* is used as a natural UV cut-off consistent with non-relativistic assumption, 
discarding fluctuation modes of order of the particle Compton wavelength [l6| . The particle then evolves in 
a new effective potential (V(r + Ar)), averaged on the fluctuation distribution, which can be written a^| 

(V(r + Ar)) = jl + ii^£)!lv 2 + . . . Jy(r) = V(r) + AV(r). 

where ellipsis dots ". . . " denote terms of order higher than first order in a. The corrective term AV(r) of 
the first order in a is precisely the term giving rise to the Lamb shift 

AE n = [ d 3 r |(r|n)| 2 AV(r) = ^-^\ g(—) (n|VV(r)|n). 

Comparison of this result with Eq. ^ shows that the IR cut-off k can be identified to (\E m — E n \) in the 
Bcthc result. This also means that k = (\E m — E a \) should not depend on the quantum numbers n (ljj . 



3. The physics of a spherical semiconductor quantum dot 

First, a QD may be considered as spherical if it possesses an aspect ratio, defined as the ratio of its longest 
over its shortest axe, smaller than 1.1. For example, according to the growth scheme used to synthesize 
CdSe nanocrystals, one can produce spherical QDs of radius in the range of a few tens A [32[ or even less 
[33l ] . as well as rod-shaped QDs, so-called nanorods, having an aspect ratio up to more than 10 [34 ] . The 
general issue of the shape of semiconducting QDs, and how QDs of different shapes should be experimcntaly 
obtained by appropriate growth technique, is extensively treated in Ref. (35| . 

Next, it is admitted that, under reasonable physical assumptions, semiconductor QDs may be described 
by EMA models in a first approximation. The most important one is the parabolicity of the semiconductor 



band structure 2^. This assumption may be efficiently amended in several ways. A review presenting 



the usual possibilities to do so is given in Ref. [111] . Here, we use the standard EMA model presented in 
Ref. 26|, as a first attempt to apprehend the Lamb effect in semiconductor QDs. Electrons and holes are 
assumed to be non-relativistic spinless particles of effective masses m* h , confined in an infinite spherical 
potential well 

if < r < R, region I ; 



V°°(r) 



00 if r > R, region II. 



The electron-hole Coulomb interaction is to be taken into account subsequently by the Ritz variational 
principle. But, they are also consequently isolated from the insulating surrounding of the QD. The electro- 
magnetic field amplitude should not then exceed some threshold, so that charge carriers would not leak out 
by tunnel effect. This working assumption is known as the weak field limit. Such model leads to an overes- 
timation of the electron-hole pair ground state energy for small QDs, which can be corrected by restoring 
a confining finite potential step of experimentally acceptable height [36|. Other QD models with parabolic 



confinement [371 l38j or parabolic potential superimposed to an infinite potential well [39j exist. But, the 



concept of QD size is no longer well-defined since their eigenfunctions are delocalizcd. 



3 Even if this Taylor expansion seems to converge only for sufficiently smooth potential V(r), it is summable in general. 
However, a fully quantum treatment, i.e. both particle and electromagnetic field dynamics are quantized, is needed. It has 
been then shown that the averaged effective potential (V(r + Ar)) can be obtained from the bare potential V(r) by applying 
the well-defined differential operator If 6l 

((Ar) 2 ) „3 

<V(r + Ar))=e 5 v V(r). 
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3.1. Interactive electron-hole pair EM A model 

The electron- hole Coulomb interaction Vc(r h) is taken into account in the total Hamiltonian 



H = H C + H h + V c (r ch ), 

where k = Aire, e denotes the semiconductor dielectric constant, r e h the electron-hole relative distance, and 
v 2 

iJ 6i h = — „ c ; h + V°°(r Cj h) the electron and hole confinement Hamiltonian. Without loss of generality, the 
semiconductor energy band gap E g may be set equal to be zero for convenience. Electron and hole, as 
decoupled particles, have wave functions of the form 

^ m (r e ,h) = R]™{r c , h )Yr{e c , hl ^) 

where I e N, n e N \ {0} and m £ [— Z,Z], Y™(0, y?) arc spherical harmonics and ]i(x) spherical Bessel 

f 1 if r € A 

functions of the first kind, xa(^) = 1 q otherwise ^ 6 rac ^ a ^ characteristic function of the set A. Finally, 

{kf%} l is the wave numbers set, defined as the n th non-zero root of ji(x), from the continuity conditions at 
r = R [23[. The respective energy eigenvalues for electron and hole are expressed in terms of {kf£}i n as 

pc.hoo _ i^ln ) 

ln " 2m* h R 2 ' 

The interplay between quantum confinement energy, scaling as oc R~ 2 , and Coulomb interaction, scaling 
as oc R~ 1 , can be described by the ratio of the QD radius R over the Bohr radius of the bulk Mott- 
Wannier exciton a* = , /z being the reduced mass of the exciton. As an exact analytic solution is beyond 
reach, two regimes are discussed in Ref. (25^ . 

3.2. Strong confinement regime 

In this regime, valid for R < 2a* , the electron- hole relative motion is so affected by the infinite potential 
well, that exciton states should be considered as uncorrelated electronic and hole states. The Coulomb 
potential is considered as a perturbation with respect to the infinite confining potential well. A variational 
approach is used to obtain the electron-hole pair ground state, with the following trial function 

0°°(r c , r h ) = ^ o (r c )V^o(rh)0rci(rch), (3) 

which is a product of non-interacting charge carriers ground state wave functions with an interaction wave 
function of the form <fi Tf} i (r e h) = c~^ Tcl ', a being a variational parameter. The electron-hole pair energy 
minimization selects the value o~o = and yields the ground state energjjf] 

E g a ** = E£-A^--4B a E*, 

where = Eq^° + Eq^° is the electron-hole pair ground state confinement energy, and E* = 2l } a * 2 ^ ne 
binding exciton Rydberg energy pjj . 



4 All constants appearing in the text and formulas are listed in appendix lAl 
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3.3. Weak confinement regime 

In this regime, valid for R > 4a*, the exciton retains its character of a quasi-particle of total mass 
AI = to* + to^. Its center-of-mass motion is confined, and should be quantized. The Coulomb interaction 
remains a perturbation to the infinite confining potential well, and the variational function </>°°(r c ,rh) 
should be kept. However, the QD size allows for a partial restoration of the long range Coulomb potential 
between the charged carriers, such that the Coulomb energy and the kinetic energy in electron-hole relative 
coordinates are of the same order of magnitude. Then, the leading contribution to the ground state energy 
of the exciton should be —E*, which may be viewed as the ground state energy of a hydrogen-like atom 
of mass fi. The total translational motion of the exciton, should be restored and contribute to the exciton 

2 

total energy by a term 2 mr? > e - the ground state energy of a free particle trapped in a region of size R. 

The exciton center-of-mass cannot reach the infinite potential well surface unless the electron-hole relative 
motion undergoes a strong deformation. The exciton should be preferably treated as a rigid sphere of radius 
7/(A)a*, where 77(A) is a phenomenologically determined function of A = ^ [25j . On this basis, to improve 
the description of the excitonic ground state, we let the center of mass motion in the variational procedure 
be represented by a plane wave </>G( r G) = e l R rTG ' ra , where tq is the center-of-mass coordinates and ctq is 
vector quantum number of unit norm |<tq| 2 = 1. The trial function </>°°(r c ,rh) takes now the form 

V oo (r c ,r h ) = oo (r o ,r h )0 G (r G ). (4) 

This leaves the exciton probability density unchanged as well as the Coulomb potential matrix element, 

2 

whereas the confinement Hamiltonian H c + mean value gets the appropriate further contribution %MW ■ 
The variational calculation yields 00 ~ 2a* _1 , and a ground state variational energy of 

7T 2 TT 2 

ch + 6/^i? 2 2M(R - v (X)a*) 2 ' 1 ' 

where terms of fourth and higher order terms in %■ are neglected. An analytical expression for 77(A) can 
now be extracted from Eq. as shown in Ref. [26|. 

3.4- Pseudo-potential-like method 

The exciton ground state energy E™£ ak of [2{|, differs from Eq. ([5]) by the term 6 * R2 , interpreted as a 
kinetic energy term in the relative coordinates. As the virial theorem in this set of coordinates should be 
satisfied, this energy should be already contained in the Rydberg energy term —E*, and therefore should 
be removed. An elegant way to do this consists in introducing a pseudo-potential 

r 2 r h S2tt 2 r 2 

W(r ch ) = W^e-^ = --^-E^o-^. 

Inspection shows that it makes contributions to the second order of the exciton total energy but not to the 

third one, while leaving the zeroth and first order terms. Higher order contributions arc interpreted as higher 

2 

order corrections to the kinetic energy of the exciton 2 mr? ■ While the amplitude W is to be fixed to get the 
correct kinetic energy ^^[^^ = — 6 * Ri j 1 + of^jj) the pseudo-potential form is not arbitrary. It is 
attractive at distances ~ a* to promote excitonic state with typical size around its Bohr radius, repulsive 
at short distances to penalize excitonic state with small size, and exponentially small for large distances in 
order not to perturb the long range Coulomb potential (26l |. 

Adding the pseudo-potential W^r^h) to the Hamiltonian Ho implies a significant decrease of the expected 
value of the exciton energy in the strong confinement regime ^ W ^^^ = — ^—CE*[l + 0(^1)}, which is 
now only valid for 2R < a* . The excitonic energy computed in presence of the pseudo-potential gets a better 
fit to experimental results in this validity domain, than those calculated without this tool. Nevertheless, the 
divergence for very small QD size still persists as a relic of the infinite potential well assumption. 
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4. Lamb shift of a particle confined in a spherical potential well 



This section deals with the Lamb shift for a particle of mass m* and of charge ±e, i.e. q = ±1, confined 
in a theoretical spherical QD with a spherical infinite potential well V°° (r The difficulty then resides in 
the determination of the Poisson equation satisfied by this singular potential. To overcome this issue, we 
will replace it by a potential step of finite constant height V > 

T/-/ \ T/ / \ f if < r < R, region I ; 

V(r) = V X]Rt00[ (r) = ^ y . f f ^ rcg . onIL 

In the limit V — > oo, the Lamb shift undergone by an energy level of the particle in the infinite potential 
well is supposed to be the finite part of the expansion in powers of V of the Lamb shift undergone by the 
level with same quantum numbers of the particle in the finite potential step, given by Eq. ([1} 



AE lnm = —log — (Vw|V "V(r)|Vw), 

OTTtTi \ K J 

where, for fixed quantum numbers I £ N, n e N \ {0} and m £ [— the wave function "0Znm(r) is the 
eigen function of the Hamiltonian Ho of the particle confined by the finite potential step V(r) instead of the 
infinite potential well V°°(rj3. 

4-1. Lamb shift 

In this new formalism, it becomes possible to find the Poisson equation satisfied by the potential step 
V(r), at least in the distributions sense. The Lamb shift of any quantum state of the confined particle is 
therefore obtained as 

Snm* 2 R 2 



AiC = -^3^1og( — )<0. (6) 



More details on this calculation is given in appendix \B.1\ As expected, because of the spherical symmetry 
of the potentials V(r) and V°°(r), the Lamb shift is independent of the azimuthal quantum number m. 
Furthermore, it is negative, since, as we shall see later, the IR cut-off satisfies the constraint n < m* . Hence, 
the Lamb effect brings down the energy levels of the particle, instead of raising them up, as in real atoms. 
This is a remarkable characteristic of the Lamb shift which, to the best of our knowledge, is found for the 
first time. 

In appendix \B.2\ we present an alternative calculation of this Lamb shift. It is based on a regularized 
but non-resumed version of Eq. ([1]), and leads to a Lamb shift, identical to Eq. Despite an identical 
regularization method, no divergent terms appear in the limit V — > oo, which is most satisfactory from 
the conceptual point of view. This corroborates the fact that such a term is not physical, when a trapped 
particle is confined by the infinite potential well V°°(r). As we will see later, this property appears to be 
fundamental. Finally, this second calculation validates the regularization of the confining potential by an 
intermediate finite potential step. 

4-2. Infrared cut-off k 

The evaluation of the IR cut-off k = (\Ef% — E°°\) is also problematic. In the Bethe approach, the IR 
cut-off explicitly depends on the considered quantum state, which is not the case in the Welton approach. 
In real atoms, it has been shown that the IR cut-off introduced by Bethe in his original article is actually 



of about the same order of magnitude independently of the energy levels |17| . Therefore, to introduce a 



5 The notation X°° = lim X is adopted, for all physical quantities. This is the reason why we will choose to note the infinite 

potential well y°°(r), and the finite potential step V(r), in sections l4l and l5l 

6 For more details, one can see appendix [B] In particular, ipi nm (r) is a linear combination of the spherical Besscl function 

jif^jf-rj and of the spherical Hankcl function h^Hi^^-rj, ki n and K\ n being the wave numbers in the regions I and II 

respectively. 
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appropriate IR cut-off k, independent from quantum numbers of the considered energy level, as required 
by the Welton approach, we suggest to define it by the formal average of the quantum state dependent IR 
cut-offs of the Bethc approach over all possible quantum numbers [171 ] 



, |poo > E mnm \E%-E%\ 1 E^(2 ? + l)(2/ + l)|(^) 2 -(fc-) 2 | 



As quantum numbers are not limited from bellow, and are infinitely many, the sum of Eq. ([7]) is infinite. 
To confer a rigorous meaning to Eq. ([7]), a rcgularization method should be prescribed. Actually, the 
number of terms in the sum is finite, since the trapped particle, being excited through its interaction with 
the electromagnetic field, can have access only to a finite number of states. As the highest accessible level is 
limited by the finite field energy En m = 2 m*R' i ; there is a way to define a UV cut-off for authorized wave 
numbers. In practice, we make use of the so-called Feynman rcgularization method, presented in appendix 
[Cl which leads to the choice 

= 7?r2 
K ~ Ylm*R? ' 

Thus, k is of the order of magnitude of the ground state energy of the confined particle, contrary to the case 



of real atoms 17J. Because of the state degeneracy, such behavior is a posteriori expected. Terms, whose 
relative importance is the largest in Eq. (O, are terms of the order of magnitude ps \{k^ +1 ) 2 — (fcj^) 2 |. The 
asymptotical behavior of the Bessel functions roots yields \{kj^ +1 ) 2 — {kf%) 2 \ « 7r 2 (2n + 1), for sufficiently 
large n € N \ {0}. By regularizing sums over n, since sums over / do not play any role, we deduce that 



n 



" m'R 2 ' 

The non-relativistic condition in the IR and UV cut-offs n(R) < m* imposes that there exists a lower 
potential radius bound i?* nin = § ~ 2, 399 X*, of the same order of magnitude of the particle reduced 

Compton wave length = m*- 1 . Thus, in a theoretical QD of radius R smaller than -R* lin , the confined 
particle should acquire at least a confinement energy E^l of the order of magnitude of its rest mass energy 

2 2 

TO * K 2m* r*' 1 — 2m* r 2 = ^oi > wn i cn w ^ explicitly contravene to the non-relativistic assumption. Finally, 
there exists another particular potential radius -R„ lax = y^eR^^ ps 3, 956 ~X*, for which the Lamb shift of any 
energy level represents a maximal fraction of it, i.e. = 0. 



dR E? 



R=R* 



5. Gauge invar iance and Lamb shift 

As an observable effect, the Lamb shift is indeed gauge invariant. Welton argument is clearly gauge 
invariant, since the electromagnetic field only intervenes through its energy density, a gauge invariant quan- 
tity. Bethe argument depends however on the gauge, since the interaction Hamiltonian H iat oc A • p in the 
Pauli-Fierz Hamiltonian Hpp, used as perturbation, is written in the Coulomb gauge. Therefore, it is useful 
to check gauge invariance in calculations of the Lamb shift of the energy levels for a theoretical QD. To this 
end, we study the Lamb effect in the gauge imposed by the electric dipole approximation. 

In this approximation, the electric field spatial variation is negligible over typical particle distances — 
i.e. E = E(t, 0) and the related scalar potential is A°(t, r) = — r ■ E(t, 0) — , and that there is no magnetic 
field B = 0, hence no vector potential A = 0. The system Hamiltonian is then of the form 

Hp F — H a + H cm + eiJj^, 

where H[' nt = —qr ■ E(i, 0) is the new interaction Hamiltonian of the particle with the electromagnetic field. 
In particular, it is obvious to observe that H[ nt and i//^ t contribute in the same way to the classical action 
of the particle-electromagnetic field system, especially when quadratic terms in A are neg lectecfl. Thus, 
quantization in the gauge fixed by the electric dipole approximation makes sense. 



7 By a simple integration by parts, one can show that J dt Hi nt =fdt H[. 
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We now go over Bethe argument with the dipole approximation. The second order correction is 



AE n = -|V^|(m|r|n)| 2 j f dk[k 2 ~ (E m - E n )k + (E, 
= -3^ 9 2 E l(m|r|n)| 2 j — -(£ m - + 

The two integrals in this expression have respectively a cubic and a logarithmic UV divergence. Thus, the 
second integral is the one which gives rise to the Lamb shift in the dipole approximation. It is natural to 
introduce m* as UV cut-off and k as IR cut-off to regularize these integrals. Rcnormalization arguments of 
the original Bethe article are used to get the corrective term 

A£ n = -!Vm* 3 (n|r 2 |n) + f q 2 m* + logf— \n|VMr)|n>. 

97T Z7T 07T TO* Z \ K J 

The third term in AE n is recognized as the Lamb shift for a state |n), as expected from the Welton approach 
or the Bethe approach in the Coulomb gauge. The second term is simply a numerical constant, which can 
be omitted using an additive rcnormalization argument. The first term is explicitly gauge dependent, since 
(n|r 2 |n) oc (n\H[^ t \n) = J2 m I ( m l^intl n ) I 2 - So, in a general way, the Lamb shift in the gauge fixed by the 
electric dipole approximation is identical to the one computed in the Coulomb gauge. 

This should be explicitly verified in the case of a particle confined in the infinite potential well V°°(r). 
As described in more details in appendix \B.3l we show that the regularized version by the finite potential 
step V(r) of Eq. (j8j), but non-resummed by the use of the closure relation, leads to the same Lamb shift than 
the one predicted by Eq. ^ for any quantum state of the confined particle. As mentioned in subsection 
\4-l\ this alternative calculation of its Lamb shift also eliminates divergent terms in the limit V — > 00. This 
reinforces once again the argument according to which the term scaling as oc VV, which appears in the 
calculation method presented in subsection \4-l\ is an artefact of the regularization method, and should be 
dropped. 



- E a fm* - (E m - E n f \og(—)\. (8) 



6. Comments on the Lamb shift of a particle confined in a spherical potential well 

The observed negative sign of the Lamb shift for a particle confined in a spherical potential found in 
subsection \4- 1\ as opposed to positive Lamb shifts known for hydrogen-like atoms, calls for a comment. We 
first present a heuristic argument, and then back it up with a mathematical approach. 

6.1. A heuristic view 

In Quantum Electrodynamics (QED), the quantized electromagnetic field is represented by an assembly 
of harmonic oscillators, which are known to have a non-zero ground state energy, often called the zero-point 
energy. Even in absence of any external radiation field, non-zero electromagnetic fields are present as arising 
from the zero-point energy of all the associated oscillators. The effect of these vacuum fluctuations on a 
charged particle is to smear out its position causing a spread of charge and mass in space, characterized by 



a typical radius of the order of magnitude of the deviation ((Ar) 2 ) 41 

In hydrogen-like atoms, the spreading effect generates a correction to the Coulomb potential. The 
electron charge is distributed over a larger volume of space and less of its charge experience the nuclear 
Coulomb potential, resulting in a reduction of its binding energy to the nucleus. Electron s-states, being 
closer to the nucleus than other states, would experience a larger reduction of the binding energy. Hence, 
the major effect of the vacuum fluctuations is to shift s-states upward in energy relative to p-states. This 
energy difference is predicted by QED to be of about +1086 MHz for the hydrogen atom 2s state, which is 
the leading contribution to the observed Lamb shift |4lj |. 

The situation is different in a theoretical QD, since no force acts on the charged particle inside the QD 
the apart infinite wall reflection at its boundary. Thus, the observed effect is not due to an electric charge 



spreading but to a mass spreading. In this situation, the total energy which was concentrated in the kinetic 
energy of the point-like particle is now transferred to the energy of a spatial mass distribution, which splits 
into center of mass motion and relative motion. On the basis of total energy conservation, one should expect 
then a reduction of the center mass motion energy, an effect which is opposite to the one observed in atoms. 

Among the secondary contributions to the Lamb shift, it is the vacuum polarization which is of largest 
order. Virtual pairs of electrons and positrons are continuously created and annihilated in vacuum. 

Inside an hydrogen-like atom, virtual electrons are attracted to the nucleus, while positrons are repelled. 
The result is that the nucleus charge is screened by the virtual electrons. But, a bound electron inside the 
cloud of virtual electrons will experience a stronger attraction to the nucleus than predicted by the standard 
Coulomb potential, and hence will be more tightly bound to the nucleus. Since s-state electrons are closer 
to the nucleus than other electrons, they will experience a larger increase in binding energy, which can be 
evaluated by QED to be of about —27 MHz for the hydrogen atom 2s state 41 1. 

But, inside a QD, the charged particle appears to move in a polarized medium in the mean. Thus, its 
motion will be kinematically enhanced, and the corresponding energy level raised a bit. In short, the Lamb 
shift in QDs is of very different nature as compared to usual Lamb shift in atoms. 



6.2. A mathematical approach 

In the calculation of the Lamb shift, the main object is the series appearing in Eq. (TTJ), which we write 
for an hydrogen-like atom as Ylijk \(^ijk\p\^nlm)\ {Ei — E„}, where 'l'nZm(r) and E n are the eigenfunction 
and the corresponding energy for a state labeled by the quantum numbers n <G N \ {0}, I € [0, n — 1] and 
m e [— I, This series is indeed convergent because of the closure relation and because E n oc -\, implying 
that E n > 0. Its sum can be shown to be ^( , I / , i im|V 2 Vc( r )| v I / nZm) oc l^nZmWI 2 > 0. The same reasoning 

noc 

holds for the series J2 ljk \{%jk\r\^ n lm)\ 2 (Ei - E n f = ^(*„Zm|V 2 Vb(r)|# niTO ) appearing in Eq. ©. 
This shows why the atomic Lamb shift is positive. 

For a QD, the two relevant series Y^ijk \(ipijk\p\ipinm)\ 2 {Eij - E hl } and Y^ijk \(^ijk\A^inm)\ 2 (Eij-Ei n ) 3 
have a completely different behavior. They are indeed not well defined, even if the regularization finite 
potential step V(r) is introduced, since the series J2ij k \('4'ijk\p\' l Pnim)\ 2 Efj (p — 1,2,3) are divergent, 
because Ei n oc kf n > oo. Yet, we have to extract a finite value from these formal series and there 

l,noo 

is no way a priori to predict the sign of the finite terrrj^|. To confer a proper meaning to these series, 
divergences should be removed by using the definition of the eigen-energies, i.e. -Solemn) = Ei„ \ipi nm ), and 
the closure relation, which stays valid because the eigenfunctions ipinmi?) are orthonormalized. This yields 
T,ijk l(V^fc|p|Vw)| 2 {Eij-Ein} = -(ipi nm \y-H a V-Ei n V 2 \^i nm ) and^-fc \{iptjk\r\^i nm )\ 2 (E tj -Ei n ) 3 = 
("i/'ir l m|[flo,r] t -[H Ql [H ,r}}\ipi nm ), as shown in appendixes EJ and EJ 

The validity of this approach is confirmed since their expansion in powers of the potential step height 
V does not contain divergent terms in the limit V oo. It also allows to understand why the expansion 
in powers of V of the matrix clement (i.'inm\^ 2 V(r)\ipi nm ) contains a divergent term behaving as oc y/V. 
Intuitively, the usual commutator relations are no longer valid, because the confined particle can have only 
access to a finite spatial region with a probability arbitrarily close to one, in the limit V — > oo. Thus, the 
effect of the surface at the boundary r = R has not been taken into account satisfactorily. If the commutation 
relations are not used, this problem vanishes. 

The sign of the Lamb shift of the confined particle in the theoretical QD is then easily predictable. 
Sections |4] and [5] suggests that it is the behavior of the radial wave function on the sphere r = R, which is 
responsible for the Lamb shift. It can be seen, from expressions of appendix [Bl that . since the contribution 
due to the Dirac distribution 5(r — R) vanishes in the limit V — > oo because of the continuity condition 
= 0, the contribution due to the Dirac distribution derivative <5'(r — R) is then the crucial one. 



8 For more details, one can refer to chapter VII in Ref. [3lll , 

9 This is a well known problem in divergent series regularization. For example, the Feynman regularization method introduced 
in appendix \B. 3\ allows to show that ^2 p>1 p 2 = — < 0. This negative extracted value is coherent with the value of the 
analytical extension of the Riemann function £ (z) = ^Z p >i to C \ {1} IS , 
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Since Rj^'(R) 7^ and ^ 0, using standard properties of the spherical Bessel functions in a 

neighborhood of one of its roots, Rj°°'(R) and R]^"(R) are of opposite sign. From this fact, it follows that 
the QD Lamb shift is negative. 



7. Lamb shift in spherical semiconductor quantum dots 

In a spherical semiconductor QD, the interaction of the electron-hole pair with the quantized electromag- 
netic field generates a Lamb shift in its energy levels. In the Welton approach, electron and hole positions 
are fluctuating. Consequently, the Lamb shift of the electron-hole pair ground state consists of the sum of 
the contributions of Lamb shift undergone by the electron and by the hole, when the pair is in its ground 
state. It will be evaluated for each confinement regime of section[3l The electron and hole respective position 
fluctuation variance 

are given by expressions analogous to Eq. ©, where k* h are electron and hole IR cut-offs. Then, in a 
semiconductor for which the fluctuations are large enough, i. e. for which the effective masses m* h are 
significantly smaller than the electron bare mass m c , there exists a possibility to detect experimentally the 
Lamb shift for a range of reasonable QD sizes. 

In this section, we perform Taylor expansions of the Lamb shift in both confinement regimes to the 
second order terms to account for the pseudo-potential W(r c h). As the contributions to the electron-hole 
pair confinement energy come only from second order terms, diagonal matrix elements of the action of the 
Laplacian on potentials are to be evaluated to this order. Therefore, the Lamb shift undergone by the electron 
or by the hole has three contributions. The first one is a boundary effect due to the confinement potential 
well y°°(r ei h), which is of the same nature as the one experienced by a confined particle in sections [4] or [5j 
The second one is due to the Coulomb interaction Vc(r e h) between the electron and the hole, analogous to 
the Lamb shift in real atoms. And, finally, the third comes from the pseudo-potential W^^). 

7.1. General considerations 

Let the non-normalized wave functions of the electron-hole pair ground state, confined by the step 
potential V(r c ^) in both confinement regimes of Eqs. ([3]) and ^ be 



(r e ,r h 



■0010 ( r e) ^oio (rh)^eh(r e h) and V( r e, rh) = i>(r e , r h )(j> G (r G ). 



Here, electron and hole ground state wave functions confined individually V'oiotTe.h) depend on the quantum 
numbers fcgi and Kq[ , subjected to k^ 2 + Kq'^ 2 = 2m* h R 2 V. The function 4>g(yg) is a pure phase factor 
and can be dropped. Since the variational parameter a appears in 0°°(r c ,rh), it should be replaced by its 
variational value (To; proper to each confinement regime. Then, the electron or the hole Lamb shift is 

AE oo = a 1 ^i^hVriV^MIf) , (^\V 2 V c (v ch )\^} , (^\V 2 W{v ch )\cj>° 



where 



3ire m* 2 h ^ K* h J\ (<^°|<£°°) 

VV(r c , h ) = V^S(r c . h -R) + 5'{r e , h - i?) j, V 2 Vfc(r ch ) - !^( 3 )(r oh ), 



and 

VMr ch )^(6-8^ + 4-t|e 

The previous mean values arc indeed confinement regime dependent, as well as the Lamb shift, because of 
the IR cut-offs k* h . Thus, we shall examine them in each confinement regime. To this end, we give some 
general results on diagonal matrix elements of the action of the Laplacian on potentials. 
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Table 1: Numerical values of material parameters in CdSo.i2'S'eo.88 ant ^ InAs, where the electron bare mass is m c 
9.11 lCT 31 kg. a a 0.529A and Ei « 13.6eV are the Bohr radius and the ionization energy of the hydrogen atom. 



Semiconductor 


in units of 


Crf5o.125eo.88 


InAs 


m* 


m c 


0.13 


0.026 


'"h 


m e 


0.46 


0.41 


e 




9.3 


14.5 


a* 


'Ki 


91.8 


595 


E* 


Ej 


1.17 10~ 3 


1.16 10" 4 



First, the case of Coulomb potential Vc(r e h) is trivially obtained since integrals over the electron (or the 
hole) coordinates are carried out over <j( 3 )(r e h)- This eliminates dependence on ctq, so that 



(<nvv c (r eh )io 



*><(r) : 



7 /d 3 r^o 4 



kR 3 



D 



is independent from the confinement regime. 

Following section [3l we have for the pseudo-potential W(r e h) 



|V 2 ^(r oh )|0 = -W^j\ 65 E - ^81 + -4^1 U / f ^^sin 2 ( 7 r 2 ;)sin 2 ( 7 rj/)sinh(I]i? a; )c- s ^, 



I? 4 



^ JoJo x V 



where £ = <r 



accounts for the exponential dependence of the pseudo-potential W^r-eh) 



The determination of the formal quantities (<f>°° |V 2 l^ oo (r ei h)|0 oc> ) is more cumbersome. As before, this 
is done using the rcgularization method of section |4l i.e. we expand the diagonal matrix element of the 
Laplacian of the hnite potential step (</>|V 2 F(r ei h)|</>) in powers of V, and assume that (a J ) oo |V 2 y oo (r Cj h)|0° o ) 
corresponds to the unique term independent from V. This calculation yields a term scaling as oc y/V, as in 
subsection^. 1 | The results of appendixes I B. 2\ and I B. 3\ suggest to drop it, because it is not physical, so that 



|V 2 V°°(r ejh )|0° 



|e- CTi? 




L 




oR 





" da; 



sin 2 (7rx) sinh(cri?a;) 



We now give the behavior of the Lamb shift of the confined electron-hole pair ground state in both confine- 
ment regimes for C(f«S , o.i2<S'eo.88 an d InAs quantum dots, with the numerical data collected in TablcQ] The 
assumption of QD sphericity turns out to be accurate for QDs synthesized in such materials. In the strong 
confinement regime, we choose to consider a QD of radius of the order of magnitude of 10A to adequately 
illustrate the properties of the predicted Lamb shift. There exists a large number of references on spherical 
CdSo.12Seo.8S or InAs nanocrystals of radius of a few tens A, which are synthesized and experimentally 
used, e.g. see respectively or 44 1. 



7.2. Strong confinement regime 

Here R < a* , Taylor expansions of the diagonal matrix elements of the action of the Laplacian on various 
potentials give the following expressions 



6°°|V 2 y°°(r ,h)|^°°) 
(<nV 2 Vb(r e h)|0°°) 



|V 2 M/(r ch )|0° 



8 p e,hoo /, p R p , R 

-&E 01 U-F- + F — 



D + F»^ + o(^ 

a* \ a* 2 



kR 3 
3 a* 2 



1 + OI 



o 



R?_ 

a* 3 
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Figure 1: Lamb shift undergone by the electron ( — ), the hole ( ) and the exciton ( — ), when the exciton occupies its ground 

state in the strong confinement regime in a. CdSo.12Seo.ss and b. InAs nanocrystals as a function of the QD radius. 




1 R (A) 



Table 2: Lamb shift undergone by the electron-hole pair ground state in Cd5o.i2Seo.88 or InAs nanocrystals a. for R = 10A or 
b. for Ft, = 30A in the strong confinement regime and c. in the weak confinement regime. 



Semiconductor 


Cd5o.i2<S'eo.88 


InAs 


E£ ( eV ) 


3.71 


15.38 




-5.52 10~ 9 


-6.1710- 7 


b- E£ (eV) 


0.412 


1.71 




-5.13 1CT 9 


-8.64 10" s 



Semiconductor 


CdSo.12Seo.SS 


InAs 


c. E* (meV) 


15.9 


1.58 


E* 


4.56 10" 7 


1.70 10" 10 



In the strong confinement regime, the electron and the hole are essentially uncorrclated, so that the IR 
cut-offs k* h should be chosen independently as, according to subsection [ 



7tt 2 



12m* R*' 



Omitting terms of the third order or higher in ^ and dropping off the exponent " 00 " , we conclude that 
the Lamb shift undergone by the ground state of the electron-hole pair in the strong confinement regime is 



where 



^strong 

77100 



16a A* 2 t 
3tt£ R 2 



log 



A£*7 = A£f°"s + AEf ons , 



R 



-F 



V.h 



R 
(f 7 



-F' 



:F" 



n *2 



<o, 



A 



and X* 



c h — iil c h auu ii mi|1 —2V 3^0 h are respectively the reduced Compton wave lengths and minimal radii 

of the electron and of the hole in the considered semiconductor. At the zeroth order in the leading term 
to the Lamb shift of the electron-hole pair ground state is negative as a symptom of the gwasi-uncorrclatcd 
confinement of the electron and the hole in the strong confinement regime. The first correction in 
particularly the contribution of the Coulomb potential Vc(r e h), is positive which recalls the behavior of the 
standard Lamb shift in real atoms. 

Figure [1] shows the behavior of the Lamb shift of the electron-hole pair ground state in spherical 
CdSo.12Seo.ss or InAs QDs. Their order of magnitude suggests a possible detection for experimentally 
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accessible QD radii. This is confirmed by tablc[2]a.-b.. since the energy orders of magnitude are equivalent 
to those theoretically predicted and experimentally observed in hydrogen atom, at least in semiconductors 
with material parameters of the same order of magnitude than those of InAs, e.g. GaAs. Here, we do 
not identify the finite potential step V(r e h) used for the rcgularization of ((jf a \V ca (Y e ' a )\(jf a ) , which is a 
calculational intermediate tool, with the real potential step, which confines the charge carriers inside the 
QD and is a physical quantity of the problem. Results of sections [4] and [5] impose that electron and hole 
masses inside or outside the QD should be identical. This means that the charge carriers are first excited, 
and then confined. However, in reality, electrons and holes are excited in a semiconducting nanocrystal. In 
this situation, electrons and holes should have different masses inside as well as outside of the QD (36^ . 

If we want to recognize the effective confining potential step at the surface of the nanocrystal with V(r e h), 
then the possibility for the electron, and to a lesser extent the hole — the notion of hole is not clearly defined 
in the surrounding insulating matrix — , to leak out of the semiconducting QD by tunnel conductivity should 
be taken into account rigourously. Since sections [4] and [5] show that the only explicit dependence on the 
step potential appears as a non-physical term in the Lamb shift expression, this refinement is superfluous, 
reflecting the coherence of the reasoning. 

Finally, in the strong confinement regime, the proposed analytical expression for the Lamb shift Ai/^amtf 
displays a linear dependence on the electron- hole pair confinement energy E^. This does not contradict 
the assumption that charge carriers are confined by a infinite potential well. Table [2] shows that its order 
of magnitude exceeds the typical height of the effective potential step at the surface of the QD. As shown 
by an analytical study of Stark effect in spherical semiconductor QDs [45| . the confinement energy itself 
should not be a relevant quantity, as far as energy shifts are concerned. Then, to validate the modeling of 
the confining potential by the infinite potential well, it is sufficient to compare the order of magnitude of 
Lamb shifts to the height of the real potential step. This is confirmed by tabic [2] a.-b., since examples of 
Lamb shifts it presents arc of the order of magnitude at most of tens /xeV, a negligible value compared to 
the effective confinement potential, which is of the order of cV f36j . 

7.3. Weak confinement regime 

For R > ira* , Taylor expansions of the Laplacian action diagonal matrix elements yield 



In this weak confinement regime, the contribution of the pseudo-potential W(r c h) to the electron-hole pair 
ground state Lamb shift cancels the second order correction of the contribution of the Coulomb potential 
Vc( r eh)- Since they both scale as <x n 2 jp = -^St, the presence of VF(r e h) allows the removal of contributions 
proportional to the electron-hole pair kinetic energy in relative coordinates, which are still contained in 
the exciton Rydbcrg energy. Finally, it is important to explain why it is not necessary to account for 
contributions coming from the confinement potential well I^°°(r h) in the weak confinement regime. Direct 
calculations show that such contributions arise as a fifth order correction term in %- to -4k 



Omitting terms of the third order or higher in and dropping off the exponent " 00 " , we conclude that 
the Lamb shift undergone by the ground state of the electron-hole pair in the weak confinement regime is 




<<nV 2 F°°(r e , h )|0 <h 





where 
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In the limit of infinite hole mass, i.e. in the limit — ^ = — > 0, in the weak confinement regime, only 

the electronic term A£"J oak contributes to the exciton ground state Lamb shift. Then, after some trivial 
manipulations, the Lamb shift of the ground state of an hydrogen-like atom of reduced mass /x in a dielectric 
medium characterized by its dielectric constant e is recognized. These observations show that electron-hole 
pai r states behave as excitonic bound states in this regime. Using known results for the hydrogen-like atom 
[171 ] , they suggest to take an acceptable approximate value for the IR cut-off as 

*s* h « 19.8£*. 

Thus, the Lamb shift undergone by the exciton ground state does not depend on the QD radius, at least up 
to the third order in This is a strong argument in favor of the validity of the pseudo-potential W{r e \ 1 ), 
which was introduced as a phenomcnological effective potential. 

Table [2] c. gives values of the exciton ground state Lamb shift in CcLSq. i2<S'eo.88 or In As nanocrystals. 
They are not experimentally accessible at present. 

7.4- Observability of the Lamb shift in spherical semiconductor QDs 

As known, the experimental observability of the Lamb shift in hydrogen atom is due to the s- and p-level 
degeneracy, when the principal quantum number is n > 2, in absence of interaction with the quantized 
electromagnetic held. The Lamb shift arises as a separation of the ns-spectral band from np-spectral band, 
while they should stay merged in absence of Lamb effect. So, how would an energy level Lamb shift be 
detected for quantum systems displaying no spectral band degeneracy such as a QD? In such systems each 
non-degenerate energy level is dressed by the quantum zero-point fluctuations of the electromagnetic held, 
forbidding the detection of the corresponding bare level. This may be the reason for the lack of articles 
addressing the Lamb effect in such structures. 

In Quantum Field Theory, the summation of the zero-point energy fluctuations yields a divergent ground 
state energy. In absence of gravity, this divergence is subtracted off in an additive renormalization scheme. 
However, a careful analysis on its volume dependence (via boundary conditions) shows the occurrence of 
a finite and observable force, known as Casimir force [46L l47|. This effect can be intuitively understood 
as follows. In vacuum, two parallel perfectly conducting squared plates of linear size L are placed at a 
separation d <C L. Since the zero-point energy fluctuations are more important outside than inside the 
plates, they are subjected to an attractive the Casimir force. 

If we adopt the Welton framework, we may view the Lamb shift as due to the particle position fluctuations 
induced by the zero-point energy fluctuations energy of the electromagnetic field. So, by placing a QD in two 
different quantized electromagnetic surroundings (in vacuum and inside a Casimir pair of conducting plates), 
one would be able to detect an energy level difference between two Lamb shifted levels. An experimental 
protocol according to which the energy levels dressed by the zero-point fluctuations energy with or without 
the two Casimir plates are to be compared. It should allow to overcome the need of degenerate energy levels, 
or of exactly computed energy levels. 

There exist some theoretical works dealing with Lamb effect of real atoms confining in a Casimir device, 
but they are based on a relativistic Bethe approach 0, [49} ■ They predict an additional shift to the standard 
Lamb shift, which depends on the separation distance between the mirrors such that it goes to zero in the 
limit d — > oo. On the other hand the Welton approach suggests, in this limit, that the Lamb shift in presence 
of the Casimir device should go to the Lamb shift in absence of such device. While the predicted additional 
shifts should be experimentally measured, they seem to not be available at present in the literature. The 
protocol suggested here is meant to reinforce the complementarity of Bethe and Welton approaches. This is 
the reason why, as a continuation of this work, we shall investigate a generalization of the Welton approach 
in order to evaluate the particle position fluctuations, in the presence of two Casimir plates. 

8. Conclusion 

Following the conventional approaches for computing the Lamb shift in real atoms, we have worked out 
an analytical expression of the Lamb shift for the ground state of a electron-hole-pair confined in a spherical 
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Table 3: Definitions, analytic expressions and approximate values of constants appearing a. in section [3] and b. in section 171 
where Oq = ago.* is the dimensionless variational parameter in the strong confinement regime. 



Name Expression Value Name Expression Value 



a. 


S 


Si(27r) ^ 


0.6720 


A 


2- 2. 

TV 


1.7861 




Bi 


2 5 

3 8ir 2 


0.6033 


B-> 


2 , 13 , S 
9 ~ l ~ 24tt 2 ~ l ~ 27T 3 " 


0.2879 




B 


B 1 + ^ 


0.6993 


B> 


AB - 1 


0.2489 




C 


i i 

3 2ir 2 


0.2827 








b. 


< 


AB' 


0.9956 


D 


2° 


1.0557 



F 1^~bXct' 0.6187 F' |_B 2 - - 0.6993 



F" BDo' 0.7350 



semiconductor QD. An explicit expression is obtained in the framework of the EMA model augmented 
by the Coulomb interaction, for both strong and weak confinement regimes. We find that, in the strong 
confinement regime and for sufficiently small but experimentally realistic QD sizes, this Lamb shift has an 
order of magnitude comparable to the Lamb shift in atoms, but is of opposite sign. 

The observability of such a Lamb shift is put to question since the energy levels are not degenerate as 
in the hydrogen atom. A way out, based on the Welton approach, is suggested. It consists in comparing 
the energy levels of the same spherical semiconductor QD in two different quantized electromagnetic field 
environments: the vacuum and a Casimir cavity. It is amusing to see that the two foremost effects which 
have served to validate QED over the years may be now reunited to demonstrate the existence of the Lamb 
shift in semiconductor QDs. 

These theoretical considerations are to be confirmed by experimental work. If the orders of magnitude 
predicted here are confirmed for the ground state of a confined electron-hole pair in the strong confinement 
regime, it would be worth while to generalize the method we have presented to any energy level of the pair. 
Therefore, the issue of how to account for the Lamb effect in modern numerical methods, such as DFT 
or ab initio methods, should be put in question. However, advances in such numerical methods are to be 
developed first. At the moment, they do not actually allow to compute energy levels of semiconducting 



nanostructurcs with sufficient accuracy [12t Il3| to observe effects of orders of magnitude of the Lamb shift. 



Finally, among several theories improving the description of semiconducting QDs by EMA models and 



their theoretical predictions [11[, let us mention a sophisticated version of the EMA, build on a reformulation 
of the multiband envelop theory (Hfl | . This is an extension of the k • p perturbation theory, which includes 
non-parabolicity in the calculation of band structure in semiconducting heterostructures |5ll . [H3 | . In the 
spherical approximation, the band structure is analytically determined in a neighborhood of a point of the 
Brillouin zone, typically its center, the so-called L-point, by second order development (5^,[H3|. It is natural 
to ask if it is possible to extend the theoretical approach presented in this paper to study the Lamb shift in 
spherical semiconducting QDs, with a standard EMA model, to the case of a multiband envelop theory, since 
all needed essential mathematical tools already exist. In particular, a orthonormal basis of eigen-functions 
of the QD Hamiltonian is explicitly known. The difficulty would therefore lie in the exact computation of 
the Laplacian matrix elements. However the feasibility of this ask does not seem, at least for the moment, 
to be attainable. 



A. Constants 

In the following tables, we sum up all appearing constants and give their approximate values. The 
function Si(x) = J Q ^ sin(i) denotes the standard sine integral. Table [3] a. (resp. Table [3] b.) gives the 
analytical expressions and approximate values of constants occurring in section [31 (resp. in section [7)1 . 
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B. Bethe approach to the Lamb shift of a particle confined in a spherical potential well 

The wave functions of the confined particle confined by the spherical finite potential step V(r) are defined, 
for quantum numbers I £ N, n £ N \ {0} and m £ ZJ, as 

Vw(r) ={Rl(r)X[o,R[( r ) +R?n(r)X[R,oo[(r)}Y l m (e,ip) 

^nj^r)x[o,i J [(r) + a„hW^^r)x [R ,oot(r-)}^ ro (^V')- 

Because the limit V — >• oo is of interest, it is possible to adjust the step potential height V, so that the 
energy eigenvalue Ei n is written as functions of the wave numbers fc;„ and Ki n in regions I and II 

1,2 K 2 

Ein = ljL ^ and V - E tn = ljL ^ > 0. 

2m*R 2 2m*R 2 ~ 

The normalization coefficients A\ n and Cz„, and the wave numbers ki„ and Ki„ are subjected to boundary 
conditions at r = R, i.e. 

f A^km) = dnh^ {iK ln ) 

\ Mnhrdl±l{hn) = iCi n Ki„\ij2i{iKin) 

and by orthonormalization conditions. In region I, the radial wave function R\ n {r) has the same form as the 
radial wave function i?J^°(r), which is then really the limit of R l ln (r), when V —} oo. The same remark should 
be made for the radial wave function R] n (r) > 0, since, by construction, K n > oo. Nevertheless, the 

Voo Voo 

introduction of the potential step allows the particle to get out the confining potential by tunnel effect, even 
if the probability for such event is exponentially small. 

B.l. Lamb shift via the Poisson Law 

The potential step V(r) satisfies the Poisson equation V 2 F(r) = V{j^5(r — R) + S'(r — i?)}, S(r) being 
the Dirac distribution and S'(r) its derivative. In the regularization method introduced in section |4l the 
unique term independent of V of the diagonal matrix element of the Laplacian of the potential step V(r) 
expansion in powers of V should be retrieved when V — > oo, since it is the only non-divergent and non- 
vanishing one. This expansion in powers of V can be explicitly done, and yields 

(Vw*|V 2 y(r)IVw) = - — \A ln \ 2 kfrji(kn)m n ) = &^—S k ?n? - ffES + 0(1). 

This suggests that the Lamb shift undergone by the particle confined by the infinite potential well V°°(r) 
is actually given by Eq. ([5]). 

B.2. Alternative route to the Lamb shift in the Coulomb gauge 

An alternative calculation of the Lamb shift undergone by a particle confined in the infinite potential 
well V°°(r) is presented here. Starting with the regularized version by the finite potential step V(r) of Eq. 
(JTJ, but for which the closure relation is not used 

&E lnm = 7^—^ W — j \{^ijk\p\^ln.m)\ {E lJ - E ln }, 

ijk 



we directly evaluate the sum J^ijk I (^'* \p\i>inm) | 2 {En - E ln } = -{ipi nm \V ■ H Q V - Ei n V 2 \%j}i nm ), and 
expand it in powers of V, as already mentioned in subsectio n^. 1\ The obtained expansion does not have 
divergent terms when V — > oo, such as terms scaling as oc \/V of appendix \B.1\ So, obtaining the Lamb 
shift for a state confined by the infinite potential well from the one confined by the finite potential step 
defined by the same quantum numbers is rigorously justified. 
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Actually, it can be shown that 

{V ■ H W - E ln } Vw(r) = R tn(R)-Jl¥n(R) L, {r _ R)+ 2 J(r _ jAynfa v)> 

from which one gets 

r>2 TP 00 

As stated, the first term in this expansion, and therefore also in the Lamb shift expansion, is the finite term, 
which leads to a Lamb shift identical to Eq. ([5]). 

B.3. Alternative route to the Lamb shift in the electric dipole approximation 

Applying the reasoning of appendix I B.2\ to the dipole approximation on Eq. (JS) , we obtain once again 
the Lamb shift given by Eq. ©. Using expressions 



AE lnm = ^ l°g(^r )J2 IOMr|Vw)| 2 (Eij - E ln f 

^ ' ijk 

^ log( — ) (V'z„ m |[ J ffo,r] t -[^o, [Ho,r]Mi nm 



37T \ K 

and computing the following commutators 

^—\d r Td r + d e r^r + d v r—^ 
m* [ r z r z sin 

and 



H ,r]ipi nm (r) = --^-i d r rd r + d e r^ + d v r " v \ij)i nm {r,9,<p), 



[H , [H 0> r]] Vw(r) = -&L^M 6 (r fl)' 
we come up with 

(Vlnmllflo.rjL [#0, [flo,r]]|Vjnm) = #T + °(!)- 

Hence, we can conclude that this Lamb shift is identical to the one given by Eq. ([5]). Let us also note that 
this method does not generate divergent terms in the limit V — > oo, as in appendix \B.2\ 

C. Feynman regularization of the infrared cut-off k 

To determine an accurate value for the IR cut-off K, the Feynman regularization method is used instead 
of the regularization method based on the existence of the UV cut-off Ki; m for authorized wave numbers 
mentioned in subsection \4-2\ Beginning with Eq. (0, a dimensionless parameter e > is introduced, and 
the regularized sum 

1 ZijinW + W + l)K fc ^) 2 - (*ff ) 2 |c-"e-*e-' e e-™ 
2m*R 2 £ 4ihl (2i + l)(2; + l)e- 8e e-J' e e-' e e-" e 

should be expanded in powers of e. The finite part of this expansion should be kept and identified with the 
IR cut-off k. Considering the Bessel functions asymptotical behavior and their roots repartition at infinity 
[40^ . for sufficiently large quantum number g£N \ {0}, but for all peff, we can assume that 

« k£ + (q-l)n » fc °? + (q- 1) = qir. 

While this approximation is not valid for any quantum numbers q € N, in the limit e — > 0, the most important 
contribution to k corresponds to high quantum numbers, i.e. it should be used for all wave numbers. Then, 

* 2 Sjnl^-jV^e-™ tt 2 f 4 3 7 ] 
K 2m*B? Y,jn e ~ 3te ~ ne 2m*R?\e 2 e 6 + ° U /' 

which yields to the definition for k, given in subsection 
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